How water droplets evaporate on a superhydrophobic substrate 
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Evaporation of water droplets on a superhydrophobic substrate, on which the contact line is 
pinned, is investigated. While previous studies mainly focused on droplets with contact angles 
smaller than 90°, here we analyze almost the full range of possible contact angles (10°-150°). The 
greater contact angles and pinned contact lines can be achieved by the use of superhydrophobic 
Carbon Nanofiber substrates. The time-evolutions of the contact angle and the droplet mass are 
examined. The experimental data is in good quantitative agreement with the model presented by 
Popov (Physical Review E 71, 2005), demonstrating that the evaporation process is quasi-static, 
diffusion-driven, and that thermal effects play no role. Furthermore, we show that the experimental 
data for the evolution of both the contact angle and the droplet mass can be collapsed onto one 
respective universal curve for all droplet sizes and initial contact angles. 



I. INTRODUCTION 

Evaporation of sessile droplets with small contact an- 
gles (< 90°) has been studied extensively. Several evap- 
oration modes have been explored: the constant contact- 
angle mode [U [2], in which the contact area of the 
droplet on the substrate vanishes; the constant contact- 
area mode [TJ [3HS]i m which the contact angle vanishes; 
and the combination of both modes [TJ [Jj [8] . A thorough 
understanding of droplet evaporation is of vital impor- 
tance for examining the drying rate [TJ [3J El |9"HTT] , the 
flow patterns observed inside drying drops [T2Tll4| . and 
the residual deposits [HI [16] . 

In early modeling of evaporating drops [31 [H [71 [T7J , the 
evaporative flux was assumed to be uniform in the radial 
direction, as it is for evaporation from a sphere. However, 
in his study of contact-line deposits, Deegan [5] argued 
that the evaporative flux from a sessile drop with a spher- 
ical cap shape is generally not uniform, but diverges near 
the edge of the drop for contact angles smaller than 90° . 
Hu and Larson [6j later used a numerical model to find an 
expression for the rate of mass loss from a drop in terms 
of its contact angle, taking this divergence into account. 
Their model applies to contact angles smaller than 90°. 

For larger contact angles, few theoretical descriptions 
exist for diffusion around a spherical-cap droplet. In [TJ, 
the rate of mass loss was expressed in terms of a series 
solution, which can be approximated in both the small 
and the large contact angle regimes. Popov [16] described 
an analytical solution for the rate of mass loss in terms 
of the contact angle, which applies to the full range of 
contact angles. However, this model has never been val- 
idated against experimental data in the large contact- 
angle regime. 

Apart from the diffusive spreading of water vapor de- 
scribed by the models mentioned above, there are other 
factors that may influence the evaporation rate (see e.g. 
[18j for an overview) . Firstly, the evaporation models dis- 
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FIG. 1. Side view (a) and top view (b) of a 8 /d droplet on a 
CNF substrate in the initial moments, (c) and (d): the same 
droplet in the last moments before being completely evapo- 
rated. Note that the contact line remains perfectly circular 
and completely pinned until almost the end of the process. 



cussed assume a stationary contact line. When the con- 
tact line is moving, dynamic effects may complicate the 
problem, for both the vapor concentration outside and 
the viscous flow inside the drop [8|. Secondly, evapora- 
tive cooling of the drop can reduce the evaporation rate 
[9TfTTj. The resulting temperature gradients on the drop 
surface can induce a Marangoni flow [T31 [2], an d can 
give rise to a Marangoni-Benard instability [19] . Finally, 
in addition to the diffusion of water vapor, free convective 
transport may play a role, increasing the evaporation rate 
[TT1 [20] , However, the influence of these factors on the 
evaporation rate still has to be confirmed experimentally. 
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FIG. 2. SEM images of the Carbon Nanofibers (CNFs) used 
as superhydrophobic substrates. Tilted side view (left) and 
augmented top view (right). 



In this paper, we describe our investigation into the 
evaporation of water droplets on Carbon Nanofiber 
(CNF) substrates; see Fig. [I] These substrates belong 
to the family of ordered carbonaceous structures: the 
graphitic planes are oriented under an angle to the cen- 
tral axis. CNF substrates can exhibit superhydrophobic- 
ity [21]. The samples used here have contact angles with 
water ranging from 150° up to 170°. On superhydropho- 
bic substrates, all evaporation modes can occur; the con- 
stant contact-angle mode is mostly observed when the 
contact angle hysteresis is low, the constant contact-area 
mode when the hysteresis is high [351 H3] • On our CNF 
substrates, the contact line remains pinned throughout 
nearly the entire experiment, hence evaporation takes 
place in the constant contact-area mode. By contrast, 
superhydrophobic substrates based on micropillar arrays 
display contact-line jumps during evaporation [251 126j . 
Since we consider pinned contact lines, we can study 
evaporating drops in almost the full range of possible con- 
tact angles (0°-150°). The rate of mass loss and contact- 
angle evolution over time are obtained experimentally for 
various drop sizes. We show that the evaporation dynam- 
ics is accurately described by the diffusion-based model of 
Popov [161, suggesting that thermal and free-convection 
effects are unimportant in our experiment. In addition, 
we show that the evolutions of the droplet mass and con- 
tact angle can be described by a universal relation, that 
is, independent of the drop size and initial contact angle. 

In Sec. [nj the experimental set-up and preparation 
of the CNF substrates are described. The experimental 
results are discussed in Sec. IIIIl The theoretical model 
for droplet evaporation adopted from Popov |16j is briefly 
described in Sec. |IV| Finally, in Sec. [Vj it is shown that 
the theoretical results are in good quantitative agreement 
with the experimental data. 



II. EXPERIMENTAL METHODS 

A. Preparation of the CNF substrates 

The droplets were left evaporating in an empty room 
[24j at a constant temperature of 23°C and a humidity 
of 30% over Carbon Nanofiber substrates; see Fig. [2] 



Carbon nanofibers were grown on oxidized silicon sub- 
strates using a Ni thin film as catalyst. A 250 nm thick 
Si02 layer was grown on p-type Si(001) via wet oxida- 
tion. On top of this oxide layer 10 nm Ta was deposited, 
followed by a 25 nm thick Ni layer. The samples were 
pretreated prior to the CNF synthesis in a quartz reac- 
tor. The substrates were placed on a flat quartz boat 
placed centrally inside a quartz reactor, and the temper- 
ature was increased at a rate of 5°C min -1 from room 
temperature up to 500°C in a N 2 (99.999%, Indugas) at- 
mosphere. During this pretreatment step, the samples 
were subjected to 20 vol.% of H2 in N2 at a total flow 
rate of 50 ml min -1 at 500 °C for 2 hours; then the tem- 
perature was increased up to 635°C. At 635 °C, 25 vol.% 
ethylene (99.95% Praxair) in N 2 was passed through the 
reactor for 1 hour, while 6.25 vol.% H 2 (99.999%, In- 
dugas) was added for the first minutes of the reaction 
time. After the reaction time, the substrates were cooled 
down in N 2 at a rate of 10° C min -1 until room temper- 
ature was reached. The CNF samples were used without 
further functionalization. 



B. Measurement of droplet evaporation 

To analyze the evaporation of droplets on CNF sub- 
strates, the droplets were observed during their total 
evaporation time and photographed at 1 s time inter- 
vals. Two synchronized cameras (Lumcncra Lml35, 1392 
x 1040 pixels) were used for this purpose, one taking side- 
view images and another one taking top- view images; see 
Fig. [T] Side-view images allowed us to compute volume 
(mass), contact angle, area, droplet radius, mass loss, 
and spreading velocity at every instant. The image anal- 
ysis was performed using a custom-made MATLAB code 
in which the detected droplet profile was fitted to an el- 
lipse. The droplets considered in this study are much 
smaller than the capillary length (which is 2.7 mm for 
a water droplet [H]), hence we can neglect flattening of 
the drops by gravity. Nevertheless, we used an elliptical 
rather than a spherical fitting. The elliptical fit allowed 
us to use three fitting parameters (two semi axes and the 
angle of the ellipse with the horizontal plane) instead of 
only one (droplet radius), thereby increasing the preci- 
sion of the determination of the volume and contact angle 
of the droplets. The ellipticity of the droplets, defined as 
the ratio between both semi axes, was always below 7%. 

The contact line of the droplets was detected automat- 
ically; the contact angles were then measured by finding 
the tangent of the ellipse at the contact line. The er- 
ror in the determination of the contact angle, based on 
the quality of the fits, was found to be below 1%. The 
volume of the droplet was obtained by calculating the el- 
lipse area above the contact line and assuming rotational 
symmetry with respect to the vertical axis, with an error 
below 10%. The rate of mass loss was computed applying 
a fourth-order finite differentiation of the ellipse volume 
over time. 
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Top-view images were used to obtain qualitative in- 
formation on the stability and circular symmetry of the 
contact line; using this information, we rejected those 
few experiments in which the contact line had a highly 
irregular shape. 

Due to the chaotic three-dimensional distribution of 
the nanofibers, the way the liquid wets the structure 
is more complex than for ordered superhydrophobic mi- 
crostructures [25] [25], for which two wetting states can be 
defined: the Cassie-Baxter state, in which the contact of 
the liquid with the substrate is minimum, and the Wenzel 
state, in which the contact is maximum. In our case, it 
is assumed that the liquid remains in a mixed state and 
that the transitions from one intermediate state to an- 
other are sufficiently smooth to be undetectable. There- 
fore, we will not use this terminology in this paper. 



III. EXPERIMENTAL RESULTS 

The droplet volume, contact angle, and radius were 
determined from the experimental data with a time reso- 
lution of 1 s. The droplet volume versus time plot, clearly 
shows nonlinear behavior; see Fig. |3ja). Hence, a model 
based on the small contact-angle approximation, which 
predicts the droplet volume to decrease linearly in time 
[5J [6], will not suffice to describe the evolution of the 
droplet volume over time. From the droplet volume mea- 
surements, the rate of mass loss of the droplet dM / dt was 
derived, as described in Sec. [n] Figure Qa) shows that 
dM/dt decreases with decreasing contact angle; hence it 
also decreases in time. Again, nonlinear behavior is ob- 
served, with a steep decline for larger contact angles, but 
this levels off to a constant rate of mass loss for contact 
angles smaller than 70°. 

During the evaporation, the contact angle of the 
droplets decreases over time from about 150° to 0°, as 
shown in Fig. [5]ja). Initially, the contact angle decreases 
slowly over time. This is followed by a more rapid, lin- 
ear decrease over time when the contact angle becomes 
smaller than approximately 70°. The initial contact an- 
gles of the droplets differ somewhat owing to irregulari- 
ties in the substrate. For comparison, not only the ex- 
perimental data, but also the predictions based on the 
Popov model are shown in Fig. [5]ja). A more detailed 
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explanation of this model is given in Sec. 

One advantage of the CNF substrates is that the con- 
tact lines of the droplets remain pinned throughout al- 
most the entire experiment. Therefore, droplet evapora- 
tion in the constant contact-area mode can be studied, in 
the absence of any contact-line dynamics. Similar behav- 
ior of the contact angle in the pinned situation has been 
reported for natural lotus leaves [23], synthetic super- 
hydrophobic surfaces with high contact angle hysteresis 
[22], and aligned carbon nanotube (CNT) samples [27] . 
Figure [6] shows that depinning only occurs during the fi- 
nal moments of the droplet life. In the depinning phase, 
the contact angle is typically smaller than 40°. Once 
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FIG. 3. (Color online) (a) Droplet volume versus time for 
initial droplet volumes of 1.6 y\ (blue filled circles), 2.1 y\ 
(red squares), 2.9 fA (green diamonds), 4.6 (A (magenta up- 
ward triangles), 6.2 pi (cyan downward triangles), and 6.9 /d 
(brown unfilled circles). The error bars are deduced from the 
elliptical fit to the data. The measurements were performed 
with a time resolution of 1 s, but for clarity we show the data 
with a 30 s resolution, (b) The dimensionless droplet mass 
plotted against the dimensionless time. The black solid line 
represents the theoretical prediction according to the Popov 
model. The experimental data is scaled according to mty. The 
time is set to at the end of the droplet life (see text). 



the droplet starts to depin, also the measurements error 
shoots up, since the contact line does not depin homoge- 
neously, and is therefore no longer exactly circular. 



IV. THEORY OF DROPLET EVAPORATION 

To describe theoretically the measured time-evolution 
of a droplct's contact angle and mass, we need to know 
the evaporative flux from the droplet surface. This flux 
depends on the rate-limiting step in the vapor trans- 
port. We assume that vapor transport by free convec- 
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FIG. 4. (Color online) (a) The rate of mass loss of the droplet 
(derived from the measured droplet volume) versus the con- 
tact angle. Colors and markers are as in Fig. [3] (b) The same 
data, but now scaled according to d6j. Predictions from the 
Popov model (black solid line) and the model of Hu & Larson 
(purple dashed line) are shown. 



tion, induced by the density difference between dry and 
humid air [20] . is negligible compared to diffusive trans- 
port. The influence of evaporative cooling of the droplet 
on the evaporation rate [TT] is also neglected. Hence, 
the vapor transport occurs mainly by diffusive spread- 
ing of the water vapor in air, and is characterized by 
diffusion time tj, = R 2 /D, with R the droplet radius 
in the plane of the substrate, and D the diffusion coef- 
ficient. The diffusion time for water vapor in air is of 
the order of 10 s. The evaporation occurs in a quasi- 
steady fashion: the timescale for diffusion is much smaller 
than the typical droplet evaporation time t e . As will be- 
come clear from the dimensional analysis presented in 
([6]), t e — p/(c s — Coo) id- In essence, t e can be estimated 
by comparing the initial droplet mass, proportional to 
the droplet density p, to the rate of mass loss, pro- 
portional to c s — Cqo, the vapor concentration difference 
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FIG. 5. (Color online) (a) The evolution of the contact angle 
over time. The experimental data (•) can be very well de- 
scribed by the theoretical model of Popov ( — ), by adjusting 
the drop radius according to its experimental value (see Sec. 
|y| . The error in the experimental data is not shown, since it 
is below 1%. (b) The same data, but with the time scaled 
according to jfjj, and set to at the end of the droplet life 
(see text). The black solid line represents the theoretical pre- 
diction according to the Popov model. Colors and markers 
are as in Fig. [3] 



between the drop surface and the surroundings. Here, 
te/td — p/(c s — Coo) is of the order of 10 5 . We do not 
take into account the Kelvin correction to the vapor pres- 
sure, since this effect is negligible for droplets of the size 
considered here. 

To determine the diffusive outflux from the drop sur- 
face, the vapor concentration field around the droplet has 
to be calculated. We follow the approach taken by Popov 
|16j . For completeness, we briefly formulate the problem 
below. 

A cylindrical coordinate system (r, z, <f>) is adopted, 
with r being the radial coordinate, z the direction nor- 
mal to the substrate, and </> the circumferential coordi- 
nate. The origin of this system is chosen such that z = 
corresponds to the substrate, and r = to the center of 
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FIG. 6. (Color online) The droplet radius versus time. Sig- 
nificant depinning of the contact line is observed during the 
final 4% of the droplet lifetime. Data are shown with 15 s 
time-resolution. During the depinning, a resolution of 5 s is 
used. Colors and markers are as in Fig. [3] 



the droplet. In this case the problem is axisymmetric, 
i.e. ^-independent. In the quasi-steady, diffusion-limited 
case the concentration field c(r, z) around the droplet is 
given by: 



0. 



(1) 



The boundary conditions imposed along the spherical- 
cap shaped droplet with arbitrary contact angle 9 are: 
(i) c = c s , the saturated vapor concentration, along the 
droplet surface; (ii) c = Coo, the ambient vapor concen- 
tration, far away from the drop; and (iii) the substrate 
is impermeable, hence dc/dz = along the substrate. 
The diffusive flux is given by J = —DVc. In our experi- 
ments, the ambient temperature was 23°C, the humidity 
H = 0.3. At this temperature, D = 24.6 x 1(T 6 m 2 /s, 
p = 997.6 kg/m 3 , and c s = 2.08 x 10~ 2 kg/m 3 (obtained 
from [55J p. 6-1, 6-191] by linear interpolation); further- 
more Coo = Hc s . 

In the limit of small contact angles, simplified solutions 
to subject to the boundary conditions (i)-(iii) can be 
used, as presented by Deegan [5] and Hu & Larson [6] . In 
our case a more advanced model is needed, since we con- 
sider droplets evaporating on a superhydrophobic sub- 
strate, with initial contact angles of approximately 150°. 
The analytical solution to the equivalent problem of find- 
ing the electric potential around a charged lens-shaped 
conductor is described in [29] ■ Popov [TS] used this re- 
sult to determine the rate of mass loss from a droplet of 
arbitrary contact angle: 



J(r) y/l + (d r h) 2 2Trrdr 



dM _ 
dt jo 

= -ttRD(c s - Coo )f{9), 



(2) 



with M the droplet mass, J the diffusive outfiux from 
the droplet surface, h(r,t) the droplet height, t the time, 



and 



/(*) = 



sin 8 



1 



cos f 



1 + cosh26»r 
sinh 2-kt 



(3) 

tanh[(?r - 6»)r]dr. 



The droplet mass can be expressed in terms of 9 by the 
geometric relation 



M{8) = pirR 



■ cos 3 9 - 3 cos 9 + 2 
3 sin 3 9 ~' 



which yields an ordinary differential equation for 
function of t 



d9 

~dt 



D(C S - Cpo) 

pR 2 



(1 + cos 8) 2 f (9). 



(4) 



as a 



(5) 



Numerical integration then gives 9 as a function of t. 
Once 9 is known, M(9) and dM/dt can be derived. 

In Fig. |3ja)-5^a) we showed the evolution of the 
droplet mass and contact angle in time for various drop 
sizes. Based on the theory described above, one would ex- 
pect a universal behavior that is independent of the drop 
size and the other problem parameters c s , H, p, and D. 
To demonstrate this, we introduce the nondimensional 
mass and time as 



M = 



M 
P~Ri 



t = 



t 



R 2 /D' 



By substituting ^ into Q-Q, we obtain 



dM 
di 

m = n'- 



de 

dt 



cos 3 9 - 3 cos 6 4 
3 sin 3 9 

(1 + cos 9) 2 f (8). 



(6) 

(7) 

(8) 
(9) 



The relations ([7|-(|9| no longer depend on the size of the 
droplets, but only on the contact angle. This implies 
that when we rescale the experimental data according to 
([6]), they should all collapse onto the theoretical curves 
described by @-([9j). 



V. COMPARISON BETWEEN THEORY AND 
EXPERIMENT 

In Sec. |IV| we explained that it should be possible 
to collapse the experimental data for all droplet sizes 
measured onto a single theoretical curve. In order to 
test this, we have to scale the experimental data accord- 
ing to ([6]). As a characteristic length scale, we would 
like to use the droplet radius. However, during the final 
moments of the droplet's lifetime, the droplet radius is 
a time-dependent quantity. Therefore, we discarded all 
data where the droplet radius was changing significantly 
(>10%) in the results that follow, and used the initial 
droplet radius for scaling. 
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The most direct prediction from the Popov model, 
which involves no time-integration, is the dependence of 
the rate of mass loss on the contact angle (|7|). Indeed, 
the scaled experimental data collapse onto a single curve, 
which is in excellent agreement with the theoretical pre- 
diction ([7]), as shown in Fig. |4jb). For comparison, the 
result obtained from applying the model of Hu and Lar- 
son [B] is also shown. Their approximation works well up 
to 9 = 90°, but for larger contact angles Popov's fully 
analytical model is required to adequately describe the 
data. 

Figure [5jb) shows that the experimental data for the 
contact angle versus (dimensionless) time follow a uni- 
versal theoretical curve for all droplet sizes measured. 
The total time it takes a droplet to evaporate depends 



applied IE1Q1]. 



VI. CONCLUSION 



on its initial contact angle, as explained in Sec. IV Since 
the initial contact angles vary somewhat, the droplet life- 
times differ. However, the experimental time is not an 
absolute measure, and we therefore have the freedom to 
set t = at whichever contact angle we want. As the 
reference point, we chose t = at the end of the evap- 
oration process, which is characterized by 9 = 0. This 
point is found by linear extrapolation from the last data 
points measured to 9 = 0. 

Once the contact angle in time is known, we can apply 
relation Q to derive the droplet mass theoretically. Ex- 
perimentally, the droplet mass is obtained independently 
of the contact angle. Therefore, the comparison between 
the theoretical predictions and the experimental data for 
the droplet mass, as in Fig. |3jb), provides a second vali- 
dation of the model. 

In the results described above, we used the experimen- 
tal data as long as the contact line remained pinned and 
hence the droplet radius remained constant. In Fig. [6j we 
showed that depinning occurs during the final moments 
of the droplet's lifetime. To construct the theoretical 
curves in Fig. [5]ja), this radius change has been taken into 
account. Time integration was performed backwards in 
time, starting from the smallest contact angle measured. 
The agreement between the model results and the ex- 
perimental data is surprisingly good, even in the regime 
where the droplet radius is changing significantly. Al- 
though the droplet radius decreases rapidly, the timescale 
over which the radius shrinks is still large -in the order 
of 100 s- compared to diffusion time (10 -2 s). There- 
fore, contact-line dynamics is still of negligible influence, 
and the quasi-steady evaporation model can indeed be 



Evaporation of water droplets on superhydrophobic 
Carbon Nanofiber substrates is studied. These substrates 
allowed us to measure the evolution of the droplet mass 
and contact angle over time, while the contact line re- 
mained pinned throughout almost the entire experiment. 
The initial contact angle was as high as 150°, and since 
it decreases to 0° during evaporation, a very large range 
of contact angles could be studied. Therefore, CNF sub- 
strates are a very useful tool to study droplet evaporation 
in the absence of contact-line dynamics. In our theoreti- 
cal analysis, we deduced universal relations for the time- 
evolutions of the droplet mass and contact angle. This 
universal scaling behavior is confirmed by our experimen- 
tal results. Since the experimental data covered almost 
the entire range of possible contact angles, we have been 
able to validate the diffusion-based analytical evapora- 
tion model presented by Popov |16) . The agreement of 
our experimental data with this theoretical model -that 
does not contain any adjustable parameters- is excellent. 
Therefore, we conclude that in our experiments the evap- 
oration is quasi-static and diffusion-driven, and thermal 
effects play no role. 

Even during the brief depinning phase, the quasi- 
steady model predicted the experimental data surpris- 
ingly well. Hence, a pinned contact line is not a strin- 
gent requirement for the applicability of the quasi-steady 
evaporation model, provided that the radius change takes 
place on a longer timescale than the diffusion. By con- 
trast, for droplets evaporating on complete wetting sub- 
strates, a quasi-static droplet profile can no longer be 
assumed and viscous effects influence the evolution of 
the contact angle over time [SJ [30! El] • It would be in- 
teresting to address intermediate cases, in which there is 
some contact-line motion, so as to establish the range of 
applicability of the quasi-steady evaporation model. 
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